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ABSTRACT 

The distributions of the stellar metallicities of K giant stars in several fields of the Galactic 
bulge, taken from the literature and probing projected Galactocentric distances of ~ 500 pc 
to ~ 3 kpc, are compared with a simple model of star formation and chemical evolution. Our 
model assumes a Schmidt law of star formation and is described by only a few parameters 
that control the infall and outflow of gas and the star formation efficiency. Exploring a large 
volume of parameter space, we find that very short infall timescales are needed ( 0.5 Gyr), 
with durations of infall and star formation greater than 1 Gyr being ruled out at the 90% 
confidence level. The metallicity distributions are compatible with an important amount of 
gas and metals being ejected in outflows, although a detailed quantification of the ejected gas 
fraction is strongly dependent on a precise determination of the absolute stellar metallities. We 
find a systematic difference between the samples of Ibata & Gilmore, at projected distances 
of 1 — 3 kpc, and the sample in Baade's window (Sadler et al.). This could be caused either 
by a true metallicity gradient in the bulge or by a systematic offset in the calibration of [Fe/H] 
between these two samples. This offset does not play an important role in the estimate of 
infall and formation timescales, which are mostly dependent on the width of the distributions. 
The recent bulge data from Zoccali et al. are also analyzed, and the subsample with subsolar 
metallicities still rules out infall timescales £ 1 Gyr at the 90% confidence level. Hence, the 
short timescales we derive based on the observed distribution of metallicities are robust and 
should be taken as stringent constraints on bulge formation models. 

Key words: stars: abundances — Galaxy: abundances — Galaxy: stellar content — galaxies: 
evolution. 



1 INTRODUCTION 

Perhaps half of the stars in the local Universe are in 
bulges/spheroids (Fukugita, Hogan & Peebles 1998) and how and 
when bulges form and evolve are crucial clues to the origins of the 
Hubble Sequence. The determination of the star formation history 
and mass assembly history of a typical bulge would provide a strin- 
gent test of galaxy formation theories. Plausible bulge formation 
scenarios range from a single high-redshift dissipative 'star-burst' 
(e.g. Elmegreen 1999), through successive merger-induced lower- 
intensity star-bursts and stellar accretion (e.g. Kauffmann 1996), 
to secular models in which predominantly stellar-dynamical effects 
at late epochs transform the central regions of thin disks to three- 
dimensional bulges via bar formation and destruction (e.g. Raha 
et al. 1991; Norman, Sellwood & Hasan 1996). The bulge of the 
Milky Way is close enough that individual stars may be stud- 
ied from the ground to within one scale-length of the centre; 
the combination of deep colour-magnitude diagrams and spectro- 
scopically determined metallicity distributions, including elemen- 
tal abundances, in principle provides both the age distribution and 
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the chemical evolution. Here we develop a model of chemical evo- 
lution of the Galactic bulge and constrain the parameter values of 
our model by comparison with the available observational data. 

There are several ways in which the bulge formation can 
be traced: a spectrophotometric analysis using multi-band broad- 
band photometry (Peletier et al. 1999; Ellis, Abraham & Dickinson 
2001) or spectral indices (Proctor, Sansom & Reid 2000) are the 
best approaches for unresolved stellar populations. However, stud- 
ies based on integrated light lack the accuracy of analyses of the 
properties of individual stars, possible for our bulge. In this cate- 
gory one should mention the study of deep stellar colour-magnitude 
diagrams (Feltzing & Gilmore 2000; Zoccali et al. 2003; van Loon 
et al. 2003) and the distribution of stellar metallicities from spectro- 
scopic surveys (Ibata & Gilmore 1995b; Sadler, Rich & Terndrup 
1996). We know from the analyses of deep colour-magnitude dia- 
grams (Ortolani et al. 1995; Feltzing & Gilmore 2000; Kuijken & 
Rich 2002; van Loon et al. 2003) that the vast bulk of the Galactic 
bulge is 'old', but the use of stellar isochrones cannot tell us how 
old to better than ~ 3 — 4 Gyr (i.e. ages in the range 10 — 14 Gyr) 
or tell us the dispersion of the distribution of ages to better than 
that. In principle, stellar elemental abundances constrain durations 
compared to lifetimes of the sources of the elements, but these are 
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difficult to obtain for large numbers of stars. We demonstrate here 
that the shape of the overall metallicity distribution provides a com- 
plementary constraint on the duration of star formation, given an 
old age. 

The available metallicity distributions for the Galactic bulge 
are reasonably well approximated by the predictions of the simple 
closed-box model (e.g. Rich 1990; Zoccalli et al. 2003). In many 
applications of this model it is a virtue that the predicted distribu- 
tion is independent of the star formation rate, but this of course pre- 
cludes conclusions about the star formation history being made by 
comparison of model predictions with observations. We have thus 
developed a model that is analytically simple, but retains explic- 
itly the timescales of gas flows and star formation. Our approach 
makes use of the metallicity distributions of K giant stars in the 
Galactic bulge. We explore a single-zone model in which no a pri- 
ori assumptions are made regarding the parameters that describe 
the infall and outflow of gas. This phenomenological "backwards" 
approach has proven to give valuable insight on the star forma- 
tion history of early-type (Ferreras & Silk 2000a, 2000b) as well 
as late- type galaxies (Ferreras & Silk 2001). Our approach comple- 
ments other studies of bulge formation, e.g. Molla, Ferrini & Gozzi 
(2000) and Matteucci, Romano & Molaro (1999), which assume a 
small set of models for a few choices of the parameters. Instead, we 
explore a large volume of parameter space, comprising thousands 
of model realizations which are later compared with the observa- 
tions. We are able to quantify how 'closed' was the proto-bulge, 
and the timescales of inflow and star formation. In §2 we describe 
our model; §3 presents the data which are used to constrain the vol- 
ume of parameter space. The comparison between model and data 
is discussed over the following two sections. Finally, §6 gives the 
conclusions to this work. 



2 THE MODEL 

The basic mechanisms describing star formation and the subse- 
quent galactic chemical enrichment can be reduced to the infall of 
primordial gas, metal-rich outflows, and a star formation prescrip- 
tion. We follow a one-zone model describing the star formation 
history in the Galactic bulge, which is a variation of the one used 
by Ferreras & Silk (2000a; 2000b) to explore the formation history 
of early-type cluster galaxies. We define our model by a set of a few 
parameters that govern the evolution of the stellar and gas content. 
A two-component system is considered, consisting of cold gas and 
stars. We adopt standard stellar lifetimes governing the ejection of 
metals but assume instantaneous mixing of the gas ejected by stars 
as well as instantaneous cooling of the hot gas component. We track 
the iron and magnesium content of the gas and the fraction that gets 
locked into stars. The basic physics related to chemical enrichment 
is reduced to: 

• Infall: The infall of primordial gas is described by a set of 
three parameters. We assume the infall rate - f(t)- to be an "asym- 
metric" Gaussian function with two different timescales (ti and Ti) 
on either side of the peak. The epoch of maximum infall is charac- 
terized by a "formation redshift" zf- Defining At = t — t(zF), we 
can write the infall rate as: 
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• Outflows: Outflows of gas and metals triggered by supernova 
explosions constitute another important factor contributing to the 
final metallicity of the bulge. We define a parameter B out which 



represents the fraction of gas and metals ejected from the galaxy. 
Even though one could estimate a (model-dependent) value of B out 
given the star formation rate and the potential well given by the total 
mass of the galaxy (e.g. Larson 1974, Arimoto & Yoshii 1987), we 
leave 73 out as a free parameter. 

• Star Formation Efficiency: Star formation is assumed to fol- 
low a variant of the Schmidt law (Schmidt 1959): 
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where p g is the gas volume density, and the parameter C e s gives 
the star formation efficiency, which - for a linear Schmidt law - is 
an inverse timescale for the processing of gas into stars. The 1.5 
exponent was chosen based on the best fit to observations from a 
local sample of normal spiral galaxies (Kennicutt 1998). This ex- 
ponent is also what one obtains for a star formation law that varies 
linearly with gas gas density and inversely with the local dynam- 
ical time, since for self-gravitating gas disks, this timescale varies 
as the inverse square root of the gas density. Changing the slope 
of this star formation dependence on gas density will alter the in- 
ferred timescales of star formation; other slopes for this correlation 
are explored in §5.3. 

The input parameters are thereby five: (ti , T2 , B out , 
C e ff , zf). The equations can be separated into one set that follows 
the mass evolution and another set that traces chemical enrichment. 
The evolution of the mass in gas and stars is given by: 
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where (f>(M) is the initial mass function (IMF). The integral E(t) 
is the gas density ejected at time t from stars which have reached 
the end of their lifetimes, tm is the lifetime of a star with mass M. 
We describe stellar lifetimes as a broken power law fit to the data 
from Tinsley (1980) and Schaller et al. (1992): 
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M t is the turnoff mass, i.e. the mass of a main sequence star which 
reaches the end of its lifetime at a time t. Finally, wm is the stellar 
remnant mass for a star with main sequence mass M: 
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The mass for white dwarf remnants was taken from Iben & Tu- 
tukov (1984). The 1.5Mq remnant mass given for the intermediate 
range is the average mass of a neutron star (e.g. Shapiro & Teukol- 
sky 1983), whereas supernovae from heavier stars might give birth 
to black holes, locking more mass into remnants (Woosley & 
Weaver 1995). The equations for the evolution of the metallicity 
of the gas and the stars are - see Ferreras & Silk (2000a; 2000b) 
for details: 

d(Z gPg )/dt =-Z g (t)il>{t) + (l-B oa t)Ez(t) (8) 
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Figure 1. Star formation history given by the model for t\ = 0.1 Gyr; 
T2 = 1 Gyr; C a g = 5; zp = 5; B ou t = 0.2. The top panel shows the star 
formation rate (if), thick line), gas infall rate (/, thin) and gas ejected from 
stars (E, dashed). The bottom panel gives the time evolution of [Mg/H] 
(dashed) and [Fe/H] (solid) in the gas, and thus in stars formed at that time. 
The inset in the top panel is the final histogram of true iron metallicities in 
long-lived stars, to be compared with the data for the bulge fields described 
in the text. 
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The yields — pm - are defined as the fraction of a star of 
mass M transformed into metals and ejected into the ISM. The 
model tracks the evolution of Mg and Fe with the yields of Thiele- 
mann, Nomoto & Hashimoto (1996) for core-collapse supernovae 
(solar metallicity progenitors) and from model W7 in Iwamoto et 
al. (1999) for type la supernovae. 

In order to estimate the rates of type la supernovae (SNela), 
we follow the prescription of Greggio & Renzini (1983), recently 
reviewed by Matteucci & Recchi (2001). We assume the progenitor 
of each SNela is a single degenerate close binary system in which 
a CO white dwarf accretes gas from the non-degenerate companion 
triggering a carbon deflagration. The rate of SNela can be written 
as a convolution of the IMF over the mass range that can generate 
such a binary system. We assume a lower mass limit of Ms m = 
3Mq in order to have a binary with a CO white dwarf which will 
reach the Chandrasekhar limit after accretion from the secondary 
star (Matteucci & Greggio 1986). Notice that this limit is somewhat 
uncertain and will be lower if He white dwarfs can give rise to 
SNela (Greggio & Renzini 1983). The upper mass limit is Mbm = 
I6M0 so that neither binary undergoes core collapse (assumed to 
happen in stars more massive than 8Mq). The rate is thus: 

Rla{t)=Aj dM<j>{M) I d///(//)V>(t-TM 2 ),(ll) 

where tm 2 is the lifetime of the nondegenerate companion, with 
mass Ma. /(m) is the fraction of binaries with a mass fraction /x = 



M2 /Mb , where M2 and Mb are the masses of the secondary star 
and the binary system, respectively. The range of integration goes 
from fi M = max(M 2 /M, 1 - 8M&/M), to (i max = 0.5. The 
analysis of Tutukov & Yungelson (1980) on a sample of about 1000 
spectroscopic binary stars suggests that mass ratios close to fi = 
1/2 are preferred, so that the normalized distribution function of 
binaries can be written: 

/( M ) = 2 1+ "*(l + 7 )^, (12) 

as suggested by Greggio & Renzini (1983), and we adopt the 
value of 7 = 2. The normalization constant A (equation II li 
is constrained by the ratio between type la and type II super- 
novae in the solar neighbourhood that best fits the observed solar 
abundances. We use the result of Nomoto, Iwamoto & Kishimoto 
(1997), namely R Ia /Rn = 0.12 to find A = 0.05, although there 
is still a rather large uncertainty in the ratio of supernova rates, so 
that Ria/Rn can be as high as 0.3 (e.g. Iwamoto et al. 1999), 
which would imply A ~ 0.12. 

We want to emphasize that the model presented here is not a 
variation of the Simple Model (e.g. Pagel 1997), as we assume that 
all of the gas in the system comes from infall as described above 
and the Simple Model usually incorporates the Instantaneous Re- 
cycling Approximation, which we do not. Our model reduces to a 
Simple Model either in the limit n , T2 — > or when we set the 
infall rate f(t) = and assume a non-zero initial gas content. 
Furthermore, the addition of SNela implies the presence of more 
high-metallicity stars (with lower [a/Fe] abundance ratios) with re- 
spect to the Simple Model, if the duration of the star formation his- 
tory is comparable to the onset time for type la supernovae. We do 
not assume any proportionality between infall or outflows and the 
star formation rate, in contrast with, e.g. Hartwick (1976) or Mould 
(1984). Our model is complementary to other work on chemical 
enrichment in the bulge (e.g. Molla, Ferrini & Gozzi 2000; Mat- 
teucci, Romano & Molaro 1999) in the sense that we let all the 
parameters controlling the SFH to vary over a wide range, only 
constraining this parameter space with the observations. Figure Q 
shows the star formation history obtained for a choice of parame- 
ters: On = 0.1 Gyr,r 2 = 1 Gyr, B out = 0.2 ,C cff = 5,z F = 2). 
The top panel shows the infall rate (/) the star formation rate (ip), 
and the gas ejected from stars (E). The lag between f(t) and ip(t) 
is caused by a combination of factors: a finite star formation effi- 
ciency, the contribution of gas from stars (E), and the power law 
used to describe the star formation rate as a function of gas density 
(equation|2j. The infall parameters t\, T2, and zf are shown in the 
figure. The bottom panel gives the evolution of [Mg/H] (dashed) 
and [Fe/H] (solid) as a function of age. One can see Mg domi- 
nates at the early stages, whereas the higher Fe yields from SNela 
make the ISM more iron rich at late times. The inset in the top 
panel is the metallicity histogram of the simulation. Infall of metal- 
poor gas provides narrower metallicity distributions by producing 
relatively more stars at later times compared with no-infall mod- 
els (such as the Simple Model). In extreme cases the infall rate 
can be balanced with the star formation rate so that the metallic- 
ity presents very narrow distributions (Lynden-Bell 1975). The low 
metallicity tail of the histogram is mainly "driven" by the early in- 
fall timescale n , which controls the buildup of metals from the pri- 
mordial infalling gas. The second infall timescale T2 has a stronger 
effect on the high metallicity part of the histogram since it only con- 
trols the infall of gas after the epoch of maximum infall, i.e. once 
the average metallicity of the ISM is rather high. The outflow pa- 
rameter Bout also plays a very important role in "modulating" the 
net yield, so that an increase in B ou t shifts the histogram towards 
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Table 1. Bulge Fields 



FIELD 


(l,b) 


N 


N sub 


M([Fe/H]) sub 


SRT96 


(+l°,-3.9°) 


268 


140 


-0.37 


IG95/#1 


(-25°, -12°) 


219 


153 


-0.57 


IG95/#3 


(-5°, -12°) 


326 


269 


-0.71 


IG95/#4 


(+5°, -12°) 


361 


271 


-0.66 


Z03 


(+0.3°, -6. 2°) 


503 


360 


-0.38 



' 1 i 1 1 1 1 i 1 1 1 1 i 1 1 1 1 r 

SRT96 




[Fe/H] 



Figure 2. Iron abundance histograms of the stars in the bulge fields used 
in this paper. Notice the offset in the distributions between the IG95 fields 
and SRT96 and Z03; a systematic offset of ~ 0.3 dex would "align" the his- 
tograms of these fields. The iron abundances are all obtained using different 
techniques, and the offsets could reflect an underlying difference in calibra- 
tions, as discussed in the text. The curve gives the expected distribution for 
the Simple Model, with a yield y = Zq . 



lower metallicities (cf. Hartwick 1976). Although only weakly, out- 
flows can also contribute to the shape of the histogram since the 
duration of the star formation stage is shortened when high out- 
flow fractions are considered, Hence, we decide to freely explore 
these three parameters (Ti,T2,B out ) and to fix the star formation 
efficiency and formation epoch to be compared for a few realiza- 
tions. We focused on three models defined by (C c s = 10,zf = 5), 
(C* c ff = 5,zf = 5), and (C c s = 10, zf = 2). The choice is jus- 
tified by the expected low fraction of young stars observed in the 
bulge, which suggests high star formation efficiencies and early for- 
mation epochs. We choose a ACDM cosmology with fl m — 0.3; 
JIa = 0.7; Ho = 65 km/s/Mpc, which implies infall maxima ocur- 
ring at around 1.2 and 3.5 Gyr for formation redshifts zf — 5 
and 2, respectively. The age of the Universe for this cosmology is 
14.5 Gyr. 



field somewhat further in projected Galactocentric distance (see Ta- 
ble 1). 

We have restricted our consideration to those fields in which 
a significantly large number of stars have derived metallicities, in- 
deed iron abundances. We should emphasize that the [Fe/H] metal- 
licities which are used to compare against our model have been ob- 
tained using different techniques. IG95 compute the 'iron' metallic- 
ity from a metallicity-sensitive Mg-spectral index (around 5200A), 
calibrated together with B — V colour (sensitive to stellar effective 
temperature) onto a [Fe/H] scale using a sample of K giant stars in 
the solar neighbourhood (Faber et al. 1985). Thus their technique 
implicitly assumes that the programme stars in the bulge and the lo- 
cal standard stars have the same value of [Mg/Fe] at a given [Fe/H]. 
As they noted, and we return to this point below, this may well not 
be the case, perhaps necessitating a new calibration of iron for their 
sample. 

On the other hand, SRT96 estimate [Fe/H] from Fe spectral 
lines plus V — I colour (again the colours are needed in order to in- 
clude sensitivity to stellar effective temperatures). The calibration 
stars are the same sample from Faber et al. (1985). Hence, a sys- 
tematic and possibly non-linear offset between these two samples 
could be expected, given the inference from high-resolution studies 
of a very limited sample, that many of the stars in Baade's window 
have enhanced magnesium, with [Mg/Fe] ~ +0.3 (Mc William & 
Rich 1994; see also Maraston et al. 2002). Finally, Z03 perform a 
metallicity estimate using only photometric data. After a careful re- 
moval of disk stars, the metallicities are computed by a comparison 
of the locus of RGB stars in the (Mk vs V — K) colour-magnitude 
diagram with analytical representations of RGB templates from 
Galactic globular clusters over a range of metallicities and thus an 
old age is assumed implicitly. As the authors caution, this anal- 
ysis may suffer from a strong systematic effect and the expected 
error bars should be larger than metallicities obtained with spec- 
troscopy. The authors further obtain an iron abundance distribution 
from their metallicity distribution by subtracting an a-element en- 
hancement of 0.2 dex for stars with [Fe/H]> —1 and 0.3 dex for 
more iron-poor stars. 



3 THE DATA 

We use the derived chemical abundance distributions from observa- 
tions of K giants in various fields towards the Galactic bulge from 
two samples: Ibata & Gilmore (1995a,b; hereafter IG95) as well as 
Sadler, Rich & Terndrup (1996; hereafter SRT96). K-giants are the 
preferred tracer since they are relatively unbiased with respect to 
age or metallicity. We also consider the recent observations of the 
bulge from Zoccali et al. (2003; hereafter Z03). IG95 target several 
fields towards the Galactic bulge along the minor and major axes, 
at projected Galactocentric distances of ~ 1 — 3 kpc, to mimic 
the long-slit spectroscopy of external bulges, whereas SRT96 con- 
centrate on Baade's window, i.e. (l,b) — (1°, —3.9°), with the Z03 



3.1 The systematics of abundance measurements 

TableQshows the fields explored in this paper along with the num- 
ber of bulge stars observed ( iV) and the number of those stars with 
subsolar metallicities (iV su b). The last column gives the median 
value of [Fe/H] for the subsolar sample. Figure [2] shows the metal- 
licity histogram of all fields considered in this paper. A significant 
offset between the IG95 fields and SRT96 can be readily seen. The 
prediction for a Simple Model of chemical enrichment with a net 
yield y = Zq is shown as a solid line; for solar elemental ratios 
this corresponds to [Fe/H]Q. The total yield is defined as: 

1 H 

V = b / dMMpu<P(M), (13) 
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where Mi is the present turnoff mass (~ IMq) and R is the re- 
turned fraction, i.e. the amount of gas returned from stars once they 
reach their endpoints, namely 

R= dM(M — wm)4>(M). (14) 

J Mi 

All histograms would peak at similar metallicities if the IG95 fields 
were shifted towards higher metallicities by ~ +0.3 dex. Notice 
that both SRT96 and Z03 have a similar value of M([Fe/H]) sub . 
Even though these two fields are not too far from each other com- 
pared to the IG95 fields, they are still separated a projected distance 
of 300 pc, i.e. around a scalelength of the bulge. The full histogram 
of SRT96 and Z03 is significantly different. 

The offset between IG95 and SRT96 can be a systematic effect 
since the estimate of [Fe/H] from Mg spectral indices and B — V 
colour uses stars with solar abundance ratios as calibrators. This is 
a perfectly valid method for local stars. However, we would expect 
stars with enhanced [Mg/Fe] - such as established for a small sam- 
ple of stars in the Galactic bulge (Mc William & Rich 1994) - to 
give systematically lower [Fe/H] if the same calibration stars are 
used. The sample of SRT96 use Fe spectral lines which should give 
a better approximation to [Fe/H] even if calibrators with solar abun- 
dance ratios are used. We can roughly estimate how much of an 
offset would be expected between these two data sets using the cor- 
rection to the metallicity suggested by Salaris, Chieffi & Straniero 
(1993). The isochrones for solar abundances corresponding to the 
corrected metallicity are equivalent to those for enhanced [Mg/Fe] 
isochrones taking the original value of the metallicity. Assuming 
an enhancement of [Mg/Fe]=+0.3 dex (Rich & McWilliam 2000), 
the correction to the metallicity obtained when the calibration stars 
have solar [Mg/Fe] would result in an offset of 0.21 dex, which is 
more or less the observed discrepancy between the median [Fe/H] 
between the IG95 fields and SRT96. 

On the other hand, this 'discrepancy' could be due to real 
astrophysics, as expected if the amount of metals driven by out- 
flows had a strong dependence on the radial position within the 
bulge. TableQshows that all three IG95 fields explored in this pa- 
per are at a Galactic latitude of —12°, whereas Baade's window is 
found at b = —3.9°, i.e. the IG95 fields are at a projected distance 
^ 1.1 kpc away from the centre, assuming a distance of 7.9 kpc 
to the Galactic centre (McNamara et al. 2000). The upshot is that 
we believe the conclusions from our analysis are very robust re- 
garding infall timescales - which pertain to relative metallicities - 
whereas the estimated fraction of gas ejected in outflows - which 
depend very sensitively on absolute metallicities - should be taken 
with care, given the possible systematic offsets between the differ- 
ent studies used. 

However, throughout this paper we shall treat all [Fe/H] esti- 
mates on an equal basis and accept them at face value. 



4 COMPARING MODEL AND DATA 

We compare the samples discussed above with our model predic- 
tion by performing a Kolmogorov-Smirnov (KS) test. The com- 
putational method is rather intensive. Furthermore, there is a de- 
generacy between infall timescale and star formation efficiency so 
that the width of the metallicity distribution depends on the ratio 
of these two parameters (see e.g. Lynden-Bell 1975). We can break 
this degeneracy by using the independent evidence of the age distri- 
bution in the Galactic bulge, specifically the low fraction of young 
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Figure 3. Probability maps at the 90% and 95% (thick) confidence levels in 
the (ti,T2) parameter space for the C c g = 10, zp = 5 (solid) and C e fj = 
10, zp = 2 (dashed) models. The dotted lines give the fraction of stars 
younger than 10 Gyr. The panels show the analysis for all the fields listed 
in table Q] The analysis is performed with the sample comprising subsolar 
metallicities. 



and intermediate-age stars inferred from deep colour-magnitude di- 
agrams (Feltzing & Gilmore 2000; Kuijken & Rich 2002; van Loon 
et al. 2003). We can thus assume a fixed (and high) value of C e ff, 
leaving the infall timescale as a free parameter. Hence, we decided 
to fix the star formation efficiency (C c g ) and formation epoch (zp) 
for three models, namely (C c ff = 10, zp = 5), (C c s = 5,zf = 5), 
and (Cog = 10, zf — 2). The remaining three parameters are ex- 
plored over a wide range: 0.01 < Ti^/Gyr < 1.5; 0.0 < B out < 
0.8.. Models with C cff = 10; n = r 2 = 0.5 Gyr form 50% of the 
stars in ~ 0.75 Gyr, compared to a longer duration of ~ 1.9 Gyr 
when the efficiency is lowered to C e g = 1. 

We computed a grid of 16 x 16 x 16 star formation histories 
and compared the simulated metallicity histograms with the data 
described above. The calibration of stars with inferred supersolar 
iron abundances is rather complicated and implies large uncertain- 
ties. Hence, two subsamples of the data were used for each field: 
one in which all observed stars were included and a second set in 
which only stars with subsolar iron abundances were included in 
the histograms. We obtained a final 3D array with each element 
representing the KS probabilities for a given choice of parameters 
(ti, T2, B ou t). Table [2] gives the 90% confidence levels for these 
parameters, using only stars with subsolar iron abundances. The 
results obtained in a comparison with the full sample is shown in 
table A KS test is an optimal statistic to be used with unbinned 
data. However, given that different statistical tests are sensitive to 
different properties of the distribution, we also performed a \ 2 test 
on binned samples. This test is rather dependent on binning. How- 
ever, the number of stars in each sample is large enough (table0 to 
make the comparison worthwhile. We computed the % 2 by taking 
10 bins in the range — 3 < [Fe/H] < +1 in the metallicity distri- 
butions of both the observed data and the models. The parameters 
we obtained using a \ 2 test were very similar to those shown in ta- 
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Table 2. Best fit parameters ([Fe/H]< 0; 90% confidence levels) 
FIELD Model n/Gyr r 2 /Gyr B out 
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+ 0.24 
-0.34 
+0.29 
-0.21 
+0.23 
-0.33 
+0.39 
-0.11 
+0.45 
-0.05 
+0.17 
-0.32 
+0.04 
-0.17 
+0.10 
-0.10 
+0.04 
-0.18 
+0.04 
-0.15 
+0.04 
-0.15 
+0.03 
-0.15 
+0.21 
-0.16 
+0.35 
-0.05 
+0.41 
-0.05 



Figure 4. Probability maps at the 90% and 95% confidence levels in the 
(ti ,-Bout ) parameter space. The same notation as in Figure|3]is used. 



estimates for the infall timescales shown in the tables do not take 
this constraint into account. However, it is worth pointing this out 
as an additional factor favouring very short infall timescales. 



bles|2|and|3|and always within the 95% confidence levels. Hence, 
we consider the best fit parameters obtained to be rather insensitive 
to the statistical estimator used. 

Figures 13141 and [5] give probability maps of the KS test at 
the 90% (thin) and 95% (thick) confidence levels. Two models are 
given: C c tt = 10, zf = 5 (solid) and C c s — 10, zf — 2 (dashed). 
Each figure shows a 2D projection of the three dimensional volume 
spanned by (ri, 72, Bout). Figure|5|shows that the model is incom- 
patible with infall timescales longer than ~ 1 Gyr. Given the high 
star formation efficiencies used in the model, this translates into a 
very short duration of the star formation stage. The results are rather 
insensitive to the formation epoch chosen. Tables|2|and|3|show that 
the best fit parameters are roughly the same for zf — 5 or zf = 2 
as long as the same star formation efficiency is chosen. This reflects 
that fact that the stellar metallicity distribution is mostly sensitive to 
the duration of star formation and does not depend much on abso- 
lute ages. A lower star formation efficiency results in longer infall 
timescales, as seen in the tables. These timescales are still shorter 
than 1 Gyr at the 90% confidence level for the more reliable subso- 
lar metallicities, except for field SRT96. 

The available colour-magnitude diagrams for lines-of-sight in 
the bulge ranging from projected distances from the Galactic Cen- 
ter of ~ 100 pc to many kpc, imply that the vast bulk of the bulge 
stars are old, with ages greater than ~ 10 Gyr (Ortolani et al. 
1995; Feltzing & Gilmore 2000; Kuijken & Rich 2002; van Loon et 
al. 2003). We incorporate these results into the comparison between 
models and data in terms of a simple fraction of stars predicted to 
be younger than 10 Gyr, shown as dotted lines in the figures. These 
mass fractions are computed for a model with C e g = 10; zf — 2 
and, where applicable, the other parameters are fixed to B out = 
(figure [3); 72 = 0.5 Gyr (figure |4j or ri = 0.5 Gyr (figure 
These parameters were chosen to represent a conservative scenario 
in the estimates of young stellar mass fractions. Hence, the figures 
show that the assumption of a stellar mass fraction no greater than 
20% in stars younger than 10 Gyr imposes a further constraint on 
the infall timescales, so that 72 ^ 0.5 Gyr would be ruled out. The 



5 DISCUSSION 

5.1 The formation history of the bulge 

Figure [6] explores the effect of varying the parameters used in 
this paper, on the stellar metallicity distribution. In all panels, the 
dashed line gives the distribution of IG95 field #1. The solid lines 
are model predictions. The top, left panel gives the best fit from 
our model to the data, corresponding to n = Ti = 0.05 Gyr; 
Bout = 0.5; Cefi = 10 (for a fixed zf = 5). The remaining three 
panels show the predicted histograms when varying any of the pa- 
rameters explored in this paper, keeping all other parameters fixed 
to the best fit values. We also show - as a comparison - the his- 
togram of the Z03 field. In the bottom, right panel, the fraction of 
gas and metals ejected in outflows is changed to B out = (keep- 
ing all other parameters unchanged). The resulting histogram cor- 
responds - to first order - to an offset towards higher metallicities, 
since more gas is allowed to be locked into subsequent generations 
of stars. The shape of the histogram will also be slightly modified 
since a low value for the outflow fraction allows a longer duration 
of star formation. On the other hand, changing the stellar yields 
would mimic a change in Bout, so that given the uncertainties in 
the Fe yields from simulations of supernova explosions (Woosley 

6 Weaver 1995; Thielemann et al. 1996; Iwamoto et al. 1999) we 
can conclude that the absolute determination of B ou t may still carry 
an important systematic offset. The bottom, left panel of figure [6] 
shows the effect of a lower star formation efficiency. The histogram 
does not change much - to be expected given the very short infall 
timescale considered - although the lower C e g tends to give lower 
metallicities. We will show below that these models can also be 
discriminated if [Mg/Fe] abundance ratios are used in the analysis. 

The shape of the metallicity distribution is strongly affected 
by a change in the infall timescale (cf. the infall solution to the lo- 
cal disk 'G-dwarf problem', Tinsley 1975). Instead of varying Ti 
and T2 separately, we show on the fop, right panel the effect of ex- 
tending the total infall timescale t/ = ri + r 2 to 1 Gyr. A more ex- 
tended infall results in a higher fraction of stars with higher [Fe/H], 
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Figure 5. Probability maps at the 90% and 95% confidence levels in the 
(t2 ,-Bout ) parameter space. The same notation as in Figure|3]is used. 



Table 3. Best fit parameters (All [Fe/H]; 90% confidence levels) 
FIELD Model rj/Gyr r 2 /Gyr B out 



IG95/#1 



SRT96 



IG95/#3 



IG95/#4 



Z03 



C10z5 
C5z5 
C10z2 
C10z5 
C5z5 
C10z2 
C10z5 
C5z5 
C10z2 
C10z5 
C5z5 
C10z2 
C10z5 
C5z5 
C10z2 



< 0.18 

< 0.31 

< 0.17 

< 0.35 

< 1.14 

< 0.28 

< 0.09 

< 0.13 

< 0.09 

< 0.09 

< 0.10 

< 0.09 

< 0.36 

< 1.15 

< 0.36 



< 0.28 

< 0.44 

< 0.27 

< 0.45 

< 0.74 

< 0.46 

< 0.10 

< 0.18 

< 0.10 

< 0.09 

< 0.13 

< 0.09 

< 0.66 

< 1.02 

< 0.61 



0.43 
0.43 
0.43 
0.11 
0.11 
0.11 
0.64 
0.64 
0.64 
0.59 
0.59 
0.59 
0.59 
0.59 
0.53 



+0.10 
0.13 
+0.11 
0.11 
+0.09 
0.15 
+0.10 
0.11 
+0.12 
0.11 
+0.08 
0.11 
+0.04 
0.11 
+0.05 
0.10 
+0.04 
0.11 
+0.03 
0.12 
+0.04 
0.12 
+0.03 
0.14 
+0.05 
0.10 
+0.05 
0.12 
+0.05 
0.05 



thereby sharpening the metallicity distribution. The prominent tail 
of the histogram observed at low metallicities shows that long in- 
fall timescales are not allowed by the observations. Quantitatively, 
figure |5| shows that star formation timescales longer than 1 Gyr 
are unlikely. Notice that field Z03 features a narrower distribution 
of metallicities, thereby favouring longer infall timescales. How- 
ever, as shown in table |2| the analysis of the (more reliable) sam- 
ple comprising stars with subsolar metallicities still discard infall 
timescales Tf ^ 1 Gyr for field Z03 at the 90% confidence level. 
It is also worth remembering that these theoretical distributions as- 
sume a very high star formation efficiency - as expected from the 
lack of young stars in the Galactic bulge. A lower value of C e g will 
imply a wider distribution of metallicities as r/ is increased. It is 
only the case of a very high C c ft along with the assumption of in- 
stantaneous mixing that gives narrow distributions when extended 
infall is assumed. 






-3 -2 -1 



HO +1 -3 

[Fe/Hl 
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Figure 6. A comparison of the effect of varying the parameters explored 
in this paper. The dashed line in all panels give the metallicity distribu- 
tion of the sample in Ibata & Gilmore field IG95#1. The solid lines are the 
model predictions. The top-left panel corresponds to t\ = T2 = 0.05 Gyr; 
Bowt = 0.5; C e ff = 10; zp = 5, which gives a good fit to the observed 
data (KS test probability 94%). The remaining three panels show the pre- 
dicted histograms when changing one of the parameters to the value shown 
keeping the other parameters fixed. In the top=right panel, Tr represents 
the overall infall timescale t\ + T2 . Notice infall timescales of 1 Gyr are 
readily ruled out because of the narrow metallicity distributions which are 
generated. The thin line in all four panels give the metallicity distribution 
of field Z03, which is narrower than IG95#1 thereby allowing for a slightly 
longer infall timescale. 



Given that the infall timescales predicted by the models are 
rather short, we decided to explore the validity of this result by per- 
forming the same test on a model which adopts the instantaneous 
recycling approximation (IRA; Tinsley 1980). In the IRA, the stel- 
lar lifetimes are assumed to be either zero or infinity, depending on 
whether the stellar mass is above or below some mass threshold, 
respectively. In that case, only stars with masses above the thresh- 
old will contribute to the enrichment of the subsequent stellar gen- 
erations. Furthermore, in this approximation we assumed type la 
supernovae do not contribute to the chemical enrichment. FigureQ 
shows the result. The probability maps of our full model are shown 
in the top panels, when comparing stars with subsolar metallici- 
ties (right), or the full sample (left) of field IG95/#1. The bottom 
panels show the result in the IRA. The infall timescales obtained 
are similar, with T\ or T2 ^ 0.5 Gyr ruled out at more than the 
90 % confidence level. Hence, the main conclusion of this paper - 
namely that star formation timescales ^ 1 Gyr are ruled out by the 
distribution of the metallicities of K giant bulge stars - is a rather 
robust statement which does not depend critically on the details of 
the adopted stellar lifetimes. 
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Figure 7. The instantaneous recycling approximation (IRA) assumes stars 
have either zero or infinite lifetimes depending on whether their masses are 
larger or smaller than some characteristic mass scale. The bottom panels 
show the probability maps when comparing the observed metallicity dis- 
tributions with a model which imposes IRA. We have also set to zero the 
contribution from type la SNe. The top panels show the results from the full 
model presented in this paper. The figure also shows the difference when 
including stars with supersolar metallicity to the analysis (left). The solid 
(dashed) contours represent the 90 (95) % confidence level. 

5.2 Linear vs Non-linear Schmidt laws 

One could expect that the short timescales obtained in the model 
presented here could vary significantly when changing the de- 
pendence between gas density and star formation. We have used 
throughout this paper a Schmidt law with an exponent of 1.5 (equa- 
tion[2j. Figure[8]shows the predicted metallicity distribution when 
the exponent is changed first to a linear, and secondly to a quadratic, 
Schmidt law. The difference between the predictions of our fiducial 
n = 1.5 model and those with n = 2 is not large. However, a lin- 
ear law (n = 1) gives a significantly sharper histogram, to be ex- 
pected since a linear law extends the period of star formation for the 
same amount of gas and star-formation efficiency. Longer star for- 
mation timescales (for fixed infall timescales) generate more stars 
with higher metallicity, thereby making the histogram narrower - 
much in the same way as more extended infall for a fixed star for- 
mation timescale. Therefore, our conclusion regarding the need for 
short star formation timescales is independent of the exponent used 
in the Schmidt law. Furthermore, a linear dependence would re- 
quire even shorter infall timescales ! 

5.3 Bulge Spectrophotometry 

We can use the star formation history that gives the best fit to the ob- 
served histograms and convolve it with simple stellar populations 
over the range of ages and metallicities predicted by the model. 
We used the latest version of the population synthesis models of 
Bruzual & Chariot (1993; priv. comm.) in order to generate an 
integrated spectral energy distribution (sed) as shown in figure [5] 
which corresponds to ri = r% = 0.2 Gyr; B out = 0.5. The sed 



has been normalized to the flux in the B-band. The resulting sed 
gives colours U — V — 1.13; V — K = 2.86. For comparison 
purposes we have considered another model with a longer infall 
timescale (n = 1 Gyr; dashed line). The colours of this model — 
U — V — 1.34; V — K — 2.99 — are redder because of the higher 
metallicities caused by longer infall timescales. The dots with error 
bars are the average and standard deviation of the U — B and B — V 
colours from the sample of 257 bulges from Sbc spirals of Gadotti 
& dos Anjos (2001). Our predicted B — V — 0.4 colour is compat- 
ible with the average and standard deviation of the observed colour 
B — V = 0.64 ±0.20 (using a subsample of bulges with nega- 
tive colour gradients). However, our predicition for U — B = 0.6 
is significantly bluer than the observed U — B — 0.19 ± 0.20. 
This may be due to uncertainties in either the chemical enrichment 
model or the population synthesis models considered. Figures |4] 
and|5|show that the value of B ou t is rather uncertain and this has a 
strong effect in the final colours. Furthermore, the discrepancy to- 
wards bluer observed colours may be caused by a systematic disk 
contamination in the observations. Even though less than 5% of the 
257 bulges observed by Gadotti & dos Anjos (2001) have colours 
U — B < 0.5, a comparison with the colours of early-type galaxies 
U — B ~ 0.6 — 0.7 (e.g. Gonzalez 1993) shows that our model 
predictions are compatible with the photometry of spheroids. 

Figure [5] illustrates the challenging task of estimating in- 
fall timescales from UBV broadband photometry alone. Only the 
spectral window around 2000A could be useful in order to rule out 
formation timescales longer than 1 Gyr. In that spectral region, the 
differences can be as large as 1 magnitude although the analysis 
would be hindered by the many uncertainties behind the model of 
star formation and chemical enrichment as well as by the uncer- 
tainties in the model predictions from population synthesis mod- 
els. The direct comparison of stellar metallicites is thereby a much 
more powerful technique to infer the star formation history of a 
bulge, stellar cluster, or galaxy. 

5.4 Simple Model and the duration of star formation 

The broad tail of the metallicity distribution at low values of [Fe/H] 
could imply that the Galactic bulge can be reasonably fit by the so- 
called Simple Model (e.g. Pagel 1997), which assumes a closed box 
system and the instantaneous recycling approximation. The beauty 
of the model lies in the ability to generate a metallicity distribution 
which is completely independent of the star formation history. The 
histogram is thereby degenerate with respect to the duration of the 
star formation process. It only depends on the total stellar yield y 
as a scale factor in the following way: 

-^^«(Z/ V )exp(-Z/»), (15) 
a log Z 

where Z is the total metallicity and M* is the stellar mass content. 
This distribution has a rather broad range of low-metallicity stars, 
which has been the reason why a Simple Model was discarded to 
explain the formation of the local Galactic disk as it generated too 
many G-dwarf stars at low metallicities (e.g. Tinsley 1980). How- 
ever, the bulge distributions shown in figure [2] are broader and so, 
more compatible with this model. 

We compared the best fits of a Simple Model with those from 
our infall model in figure [T0l We plot the Kolmo gorov-Smirnov 
probability when comparing fields IG95#1 (top) and SRT96 (bot- 
tom) with our model for a fixed formation redshift: zf = 5, 
Ti = 0.1 Gyr. The outflow parameter was chosen to maximise 
the probability for each field. We take B out = 0.5 for IG95#1 and 
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Figure 8. Metallicity distribution predicted for three different exponents of 
the Schmidt law relating the star formation rate with the gas density. The 
model corresponds to n = T2 = 0.1 Gyr; B out = 0.5; C c g = 10. 
Notice that a linear law gives a much sharper histogram since star formation 
takes longer for the same amount of gas, compared to a non-linear law with 
n > 1. 

Bout = for SRT96. Three different star formation efficiencies 
are considered, namely C c s = {5, 10, 20}. The KS probability is 
shown as a function of T2 (left) or AtsF (right), where the latter 
is defined as the time lapse during which 75% of the total stellar 
mass content in the bulge is generated. Only sub-solar metallici- 
ties are considered in the test. The horizontal line gives the highest 
probability for a Simple Model, varying the yield (y) in order to 
maximize the probability. Hence, one can see that our models give 
better fits than a Simple Model, so that a detailed analysis of the 
metallicity distribution in the bulge can be used to set limits on 
the duration of star formation. All models shown in the figure give 
star formation timescales shorter than ~ 1.5 Gyr. Furthermore, the 
[Mg/Fe] enhancement observed in bulge stars (Rich & Mc William 
2000) speaks in favour of models with a high C e g so that star for- 
mation timescales AtsF ~ 0.5 Gyr are expected. 



6 CONCLUSIONS 

We have explored a simple model describing the formation and evo- 
lution of the stellar populations of the Galactic bulge using a set of 
a few parameters. A comparison of the model predictions with the 
observed [Fe/H] distribution of K giants in various fields towards 
the bulge requires relatively short infall timescales ( £ 0.5 Gyr) re- 
gardless of the field considered (Figure[3j- Our model compares the 
resulting star formation history and subsequent enrichment with the 
distribution of stellar metallicities. Within the model assumptions 
and other uncertainties, one can relate the infall timescales with 
the duration of star formation. Hence, formation timescales longer 
than Si 1 Gyr are ruled out at more than the 90 % confidence level 
regardless of the field, statistical test, or on whether stars with su- 
persolar metallicities are excluded from the analysis (Figure^. 




- T 1= lGyr 

2000 3000 4000 5000 6000 

\(A) 

Figure 9. Spectral energy distribution (sed) obtained by convolving the star 
formation history for the best fit to Ibata & Gilmore field IG95#1 with the 
stellar population synthesis models of Bruzual & Chariot (in preparation). 
The model corresponds to n = T2 = 0.05 Gyr; B ou t = 0.5. The vertical 
lines give the central positions of the U, B, and V passbands. A second 
sed with a longer timescale (ri = 1 Gyr) is also shown (dashed line). The 
larger population of stars with a higher metallicity in the t\ = 1 Gyr model 
is responsible for the redder colours. The dots correspond to the average 
and standard deviation of the Sbc bulge sample of Gadotti & dos Anjos 
(2001). A comparison of the solid and dashed lines shows that it is very 
difficult to determine SFH parameters from integrated spectroscopy alone. 
The distribution of stellar metallicities (as presented in this paper) is a much 
more sensitive discriminator. 



Outflows during the formation of the bulge can also be esti- 
mated, although this result is very strongly dependent on the (un- 
certain) absolute stellar yields of Fe from type II supernovae and 
requires a precise calibration of the absolute metallicities of the 
observed stars, including the effect of non-solar abundance ratios. 
Taking the derived distributions as face value, even if we decided 
to track the total metallicity, Z, instead of the Fe content, the un- 
certainties in the stellar yields from intermediate mass stars are still 
large enough so that the model estimates of the fraction of gas and 
metals ejected in outflows are more qualitative than quantitative. 
Nevertheless, we find outflows may be very significant in most of 
the fields studied. Baade's window (SRT96) is the field in which 
high outflow fractions are ruled out. On the other hand, the fields 
from IG95 favour a large amount of gas ejected in outflows. In these 
fields, Bout ~ 0.5 — 0.7 give the best fit, and B ou t = seems to 
be a very unlikely scenario even if the real stellar yields or the IMF 
are far from those adopted in this paper. This difference in outflow 
fraction is in the sense expected if outflow is inhibited in the deep- 
est part of the potential well. Large outflow fractions are expected 
and predicted in some models that include estimates of the dynam- 
ical effects of feedback from stars (e.g. Arimoto & Yoshii 1987). It 
is worth noticing that these values of B ou t imply a very significant 
amount of metals contributing to the enrichment of the IGM (Ren- 
zini 1997), provided the gas escapes the overall Galaxy, rather than 
for example enriching the disk. 
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Figure 10. Comparison of the Kolmogorov-Smirnov probability as a func- 
tion of T2 (left) or Atgp (right). The latter gives the time for the bulge 
to generate 75% of the total stellar mass content. Fields IG95#1 (top) and 
SRT96 (bottom) are tested with a fixed value of B out = 0.5 (IG95#1) and 
B out = (SRT96) and ti =0.1 Gyr, z F = 5 for both fields. Three star 
formation efficiencies are chosen, namely C ff = 5 (dashed), 10 (dotted), 
and 20 (solid). The horizontal line gives the best fit for a Simple Model. 



The short timescales predicted by the analysis of the metal- 
licity distribution of bulge stars has a direct consequence on abun- 
dance ratios such as [Mg/Fe]. This ratio is a reasonably robust in- 
dicator of the duration of star formation. Solar values are achieved 
when extended star formation takes place, so that the debris from 
SNela can be incorporated into subsequent generations of stars. 
On the other hand, short-lived bursts such as the one we predict 
in this paper, translates into enhanced [Mg/Fe] in most bulge stars. 
Figure \H] shows this point. Analogously to figure [6] we give the 
model prediction for various choices of the parameters including 
the set that gives the best fit. In this case the histogram of [Mg/Fe] 
is shown. In most cases, the distribution is very similar, peaked 
at [Mg/Fe] ~ +0.25 with an extended tail which dies off at so- 
lar abundance ratios. Only the model with low star formation effi- 
ciency (bottom left) gives a significantly different histogram. A low 
Ceft implies a more extended period of star formation, generating a 
broader histogram as more stars become polluted by SNela ejecta. 
The abundance ratios observed by Rich & McWilliam (2000) on 
a sample of bulge giants using Keck/HIRES rule out models with 
low star formation efficiencies. 

The main outcome of this paper is that infall timescales 
r ^ 1 Gyr are ruled out by the observed metallicity distribution 
of the Galactic bulge. Short star formation times are also to be ex- 
pected given the observed broad distribution of the metallicities of 
bulge K giants. The preferred value stays around infall timescales 
Tf ~ 0.1 — 0.2 Gyr, which would correspond to star forma- 
tion timescales tsf ~ 1 Gyr. However, the differences between the 
iron-abundance distributions from the different fields explored in 
this paper illustrate the need for a uniform large-scale (IR) spec- 
troscopic survey, with radial velocity AND STAR COUNTS used 



Figure 11. Distribution of [Mg/Fe] abundances predicted by the models 
shown in figure |6| The histogram corresponding to the best fit is shown 
in the top-left panel as well as in the three remaining panels as a dashed 
line. Notice that a low star formation efficiency - which generates a similar 
[Fe/H] as the best fit as seen in figure^]- generates a very wide [Mg/Fe] 
distribution. 



as an aid to the statistical separation of bulge and foreground disk 
AND STELLAR HALO. 

A robust calibration of metallicities - which should account 
for variations in the [a/Fe] abundance ratios - is crucial in the quan- 
tification of the amount of gas ejected in outflows during the for- 
mation of the bulge. Furthermore, an accurate distribution of abun- 
dance ratios such as [Mg/Fe] (figure ITTI would pose strong con- 
straints on the duration of the star formation burst which gave way 
to the Galactic bulge. 

Our technique is quite robust and complementary to a com- 
parison of the stellar colour-magnitude diagrams of old popula- 
tions. The latter can only constrain the star formation timescale 
to a few Gyr due to the crowding of the isochrones at old ages. 
The model presented here sets strong constraints for bulge forma- 
tion in semi-analytical models, the latest of which predict forma- 
tion timescales for bulges significantly longer than these timescales 
(Abadi et al. 2002). Furthermore, secular evolution models predict 
formation times which extend over many dynamical timescales. 
Hence, formation scenarios over times ~ 0.1 — 0.2 Gyr, as pre- 
sented here, would be in conflict with bulge formation from bar in- 
stabilities. Our results for the buildup of the Galactic bulge favour 
formation scenarios in which a strong starburst quickly converts all 
gas available into stars in a few dynamical timescales (Elmegreen 
1999). 
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